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A theory describing a one-dimensional Luttinger liquid in contact with superconductor is developed. 
Boundary conditions for the fermion fields describing Andreev reflection at the contacts are derived 
and used to construct a bosonic representation of the fermions. The Josephson current through a 
superconductor/Luttinger liquid/superconductor junction is considered for both perfectly and poorly 
transmitting interfaces. In the former case, the Josephson current at low temperatures is found to be 
essentially unaffected by electron-electron interactions. In the latter case, significant renormalization 
of the Josephson current occurs. The profile of the (induced) condensate wavefunction in a semi- 
infinite Luttinger liquid in contact with a superconductor is shown to decay as a power-law, the 
exponent depending on the sign and strength of the interactions. In the case of repulsive (attractive) 
interactions the decay is faster (slower) than in their absence. An equivalent method of calculating 
the Josephson current through a Luttinger liquid, which employs the bosonization of the system as 
a whole (i.e., superconductor, as well as Luttinger liquid) is developed and shown to give the results 
equivalent to those obtained via boundary conditions describing Andreev reflection. 

PACS numbers: 74.80.Dm, 74.80.Fp, 73.20.Dx 



I. INTRODUCTION 

When a normal metal (N) is put in good electrical con- 
tact with a superconductor (S), superconducting order is 
induced in the normal metal over distances far greater 
than any microscopic lengthscale, either of the normal 
metal or the superconductor. This induced order leads to 
a number of remarkable phenomena, such as the Joseph- 
son effect in SNS junctions |jj and the induced Meissner 
effect in SN bilayers Q], collectively known as "proxim- 
ity effects" Until very recently, all work on such 
effects, both experimental and theoretical, has concen- 
trated on systems in which N is in the Fermi-liquid (FL) 
state. It has long been appreciated theoretically that, in 
contrast with their higher-dimensional analogues, (effec- 
tively) one-dimensional systems of interacting electrons 
are not Fermi liquids. Instead, they exhibit a number of 
possible regimes (3), among which the Luttinger liquid 
(LL) provides a one-dimensional (ID) metallic counter- 
part to the (higher-dimensional) FL state, albeit differing 
in several important respects, most notably in the ab- 
sence of single-particle excitations in the low-energy part 
of the spectrum. The basic features of LLs have been 
understood mainly in the context of ID organic charge- 
transfer and mixed- valence conductors Q). In addition, 
the prediction of the suppression of the tunneling conduc- 
tance of LLs [||J6) has stimulated the experimental search 
for Luttinger liquids in mesoscopic systems, in particular, 



in the edge-channels of fractional quantum Hall systems 
and in semiconductor quantum wires Q . 
The purpose of this paper is to address the is- 
sue of proximity effects at Luttinger- liquid/supercon- 
ductor interfaces, including the Josephson effect in 
superconductor/Luttinger- liquid/superconductor (SLS) 
junction. Our motivation is twofold. First, experimen- 
tally, such a study is relevant in view of the rapid progress 
in the fabrication of superconductor/semiconductor in- 
terfaces B, especially those with high interface-trans- 
parency (such as, e.g., the Nb/InAs interface), and also 
in view of the recently reported observations of LL-like 
behavior in GaAs quantum wires ||. Thus, the fabri- 
cation and investigation of SLS systems may reasonably 
be anticipated in the near future. Second, theoretically, 
we aim to understand the interplay between electron- 
electron interactions and induced superconducting order 
in ID electronic systems. Furthermore, one of the pos- 
sible scenarios of high-temperature superconductivity in 
oxide materials is built on the assumption of the LL-like 
character of the normal electronic state in these materi- 
als jll| . The existence of LLs in dimensions higher than 
one, however, is not yet established, in contrast to the 
ID case. Thus, a ID LL in which the superconductivity 
is induced via the proximity effect may provide a model 
system for superconductivity in 2D [|l2| . 

Our main results can be formulated as follows: (i) At 
low temperatures, the Josephson current through an SLS 
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junction having perfectly transmitting interfaces has the 
same phase- and length-dependences as in the noninter- 
acting case, the only difference being a renormalization 
of the effective Fermi velocity. The reason for this is that, 
using the bosonic language, the non-dissipative (topolog- 
ical) currents, including the Josephson current, are car- 
ried in LLs by the topological modes of the boson fields, 
which are not sensitive to the interactions. At temper- 
atures above a certain crossover value, interactions lead 
to the additional suppression of the Josephson current, 
(ii) The (induced) superconducting condensate wavefunc- 
tion in a LL in good electrical contact with a supercon- 
ductor decays as £~ 7 with the distance x from the LL/S 
interface, with 7 depending on the strength of the in- 
teractions. (7 = 1 corresponds to the noninteracting 
case, whereas 7 > 1 for repulsive and 7 < 1 for attrac- 
tive interactions.) (iii) For the case of imperfectly trans- 
mitting interfaces, the renormalization of the interface- 
transmission coefficients (via a mechanism known from 
studies of LLs coupled by weak links ||[§) ^ s renec ted in 
the renormalization of the Josephson current, which gets 
strongly suppressed in the case of the repulsive interac- 
tions. Along the way, we have also: (iv) derived effec- 
tive boundary conditions describing Andreev reflection 
at the interface with a superconductor, which we have 
then used as boundary conditions in the bosonization 
procedure; (v) determined the structure of the topologi- 
cal (Haldane) excitations in an SLS system; and (vi) con- 
firmed result (i) via an alternative approach, in which 
bosonization is applied to both the superconducting and 
normal parts of the system. 

The issue of the Josephson current through a LL 
has also been studied in a recent paper by Fazio, 
Hekking, and Odintsov (ll| for the case of poor interface- 
transmittance (see also Rcf . |T(| ) . By using the tunneling 
Hamiltonian method, it was found that the Josephson 
current through an SLS junction is suppressed compared 
to the noninteracting case. The present paper takes a 
different approach. This approach originates from work 
on SNS junctions with perfect interface-transmittance 
fTi] |l7f , in which the Josephson current was related to the 
spectrum of electronic states confined to the N region by 
Andreev reflection. Our results concerning poorly trans- 
mitting interfaces agree with those of Ref. |l| (up to 
non-universal numerical prefactors, which we have not 
attempted to calculate). 

The present paper is organized as follows. In Sec. |l| 
we derive the boundary conditions for the fermion field 
operators at the NS interface in the absence of interac- 
tions. In Sec. [II we develop a bosonization procedure for 
interacting fermions confined to the normal ID region of 
an SLS system, which makes use of the boundary con- 
ditions derived in Sec. [ij. We calculate the Josephson 
current through an SLS junction in Sec. IV. In Sec. [v|, 
we analyze the profile of the condensate amplitude in 
a semi-infinite LL connected to a superconductor. Up 
to this stage, we will have applied bosonization only to 
the normal part of the system, the presence of a super- 



conductor being implemented as a boundary condition. 
An alternative approach, in which both the normal and 
the superconductin g p arts of the system are bosonized, 



is presented in Sec. VI 



II. ANDREEV BOUNDARY CONDITIONS 

A. Andreev reflection at the NS interface: 
qualitative picture 

Electronic excitations in a normal metal having ener- 
gies smaller than the superconducting energy gap Ao suf- 
fer Andreev reflection at the interface, i.e., electron- like 
excitations are reflected as hole-like excitations, with a 
Cooper pair being injected into the superconductor, and 
vice versa Jl|,[l!| • The single-particle excitations in S are 
mixtures of electron- and hole- like states with weights de- 
termined by the self-consistency condition. In the bulk 
of N, the electron- and hole-like states are uncorrelated. 
Near the boundary, however, Andreev reflection mixes 
the electron- and hole-like states precisely in the same 
proportion as they are mixed in S, which leads to the for- 
mation of a condensate, the amplitude of which decays 
into the bulk of N. The decay-length of this condensate 
is the length Lt over which superconducting correlations 
in the motion of bulk normal-state electrons exist. In 
the case of perfect metals, Lt = Hvp/T, where is 
the Fermi velocity and T is the temperature (we choose 
units in which — 1). (The same length determines the 
thermal disruption (in the absence of inelastic processes) 
of mesoscopic phase- coherence, as is manifested in the 
phenomenon of universal conductance fluctuations |po| ) . 
For T -C A ~ T c (where T c is the critical temperature 
of S), Lt ^> £s (where £s — Tivf/Ao is the coherence 
length of S). In order to describe the influence of the 
superconductors on the N region, we now derive effec- 
tive boundary conditions that account for the Andreev 
reflection suffered by the low-energy components of the 
fermion fields at th e NS interfaces. Our strategy is as fol- 
lows: in Sec. [IB, we derive these boundary conditions 
for the case of the non-interacting electron gas in N; then, 
in Sec. Ill, we implement these boundary conditions into 



the bosonization scheme for interacting electrons. 



B. Derivation of Andreev boundary conditions 

We consider a one-dimensional electronic conductor 
(i.e., a quantum wire) of length L, adiabatically con- 
nected to superconducting leads (see Fig. [Ja). We begin 
by analyzing the ideal case, in which the single-electron 
parameters (Fermi velocities, effective masses, etc.) are 
the same in the N and S parts of the structure, the only 
difference between N and S being the presence of a pairing 
potential in S. We adopt the conventional model jf4| |T^] 
in which the pairing potentials in the leads are assumed 
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to be unaffected by the presence of N. Although this is 
a non-self-consistent approximation, it is known to re- 
flect correctly the aspects of the problem relevant for the 
present treatment [^4|-^7|. 

We temporarily replace the real 3D superconducting 
leads by effective ID leads. The profile of the pair po- 
tential is then given by (cf. Fig. 



A(x) = 




for x < 0; 
for < x < L; 
for x > L. 



(1) 



In the Andreev (semi-classical) approximation |fl8| , 
which is valid for Ao « e hkpv-p/2, the spinor of 
Bogoliubov amplitudes 



(2) 



satisfying the Bogoliubov-de Gennes equations {22J is de- 
composed into left- and right-moving components, 



w = e lkFX vf + 



(3) 



where k-p is the Fermi wavevector. The components 
w± now satisfy the (formally) relativistic (first-order) 
Bogoliubov-de Gennes equations: H^w± = ew±, with 
the Hamiltonians 



^fihvpdx A(ai) 
A*(ar) ±ihvpd x 



The full solution of these equations is obtained juj by 
finding the solutions in the N and S regions and then 
matching them at the interfaces. (In the semi-classical 
approximation, only the wavefunctions need be continu- 
ous.) The solution in N for e < Ao can be written as 



w± = 



where 



11(e) 



and 7/(e) = cos 



(4) 



(5) 



in which 72 is the Andreev reflection coefficient, whose 
phase is 77. The quasiparticle momentum tik = e/v-p sat- 
isfies the quantization condition 



72(e) 



2+i( X l-X2)„2ikd 



1. 



(6) 



where ± corresponds to two sets of energy levels O] . In 
Eq. (Q), A± are overall normalizations which, without the 
loss of generality, can be chosen to be real. Evaluating 
Eq. (Q) at x = and x = L and using Eq. @, one can 
see that at the boundaries the left and right components 
of the Bogoliubov amplitudes satisfy 



v± 
v± 



7^Ti e -*xi u±j 
72 



->X2 



u±, 



for x 
for x 



0, 
L. 



(7) 



Equations (^) describe the essence of Andreev reflection: 
the electron- like excitations (u±) are converted into hole- 
like excitations (v±), at the same time acquiring the 
phase of the order parameter (xi.,2) together with the 
phase-shift of the Andreev reflection coefficient (77) . 

In the limit e <C Ao, the phase shift r\ — * 7r/2 |2l]], and 
the boundary conditions (Q) become energy-independent. 
This enables one to derive from Eqs. (|7|) the boundary 
conditions for the real-space fermion operators tl> s (x), 
where s — | , | denotes the spin projection. These field 
operators are related to the (it, v) amplitudes via the Bo- 
goliubov transformation p2[ 



(8) 



where c s (c\) is the fermion annihilation (creation) oper- 
ator, the sum runs over all single-particle quantum num- 
bers, and the variable s takes on the values +1(— 1) for 
the f (|) spin-projections. We decompose ip s (x) into the 
left- and right-movers: 



ik-px 



(9) 



Substituting decompositions (||) and (||) into Eq. (j|), we 
obtain the Bogoliubov transformation for ip±,s- 



(10) 



The boundary conditions for the (Pauli) spinors ip±, 9 
then follow upon substitution of Eqs. ^ into Eq. (|lT 
and using rj = tt/2. After some algebra, we obtain 



,l\x=0,L ' 



or, more compactly, 



■hie 



<Xl,2,/,t 



x=0,L 



x=Q.L 



(11a) 
(lib) 



(12) 



Here, T — gC is the time-reversal operator [ p2| , with 

g = ia y = f _°! q ) ( 13 ) 

and C being the Hermitian conjugation operator. The 
presence of T in Eq. ( |12[) signals an important property 
of Andreev reflection |18| : a reflected excitation is the 
time-reversed version of an incident one. 

Further insight into the meaning of the boundary con- 
ditions ([12]) can be obtained by employing the chiral sym- 
metry of left-right fermion fields ">P±.s- In what follows, 
we adopt the methods of Refs. , in which the chiral 

symmetry of boson fields satisfying Dirichlet or Neumann 
boundary conditions was used to derive effective periodic 
boundary conditions. The right (left) field describes the 
propagation of the (formally) relativistic fermions to the 
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right (left) with the Fermi velocity. Consequently, in the 
Heisenberg representation, the space-time dependence of 
these fields is given by 



ip±, a (x, t) = ip±, s (x =F v F t). 



(14) 



Using the boundary conditions ( lla| , llb| ), one sees that 
at any instant of time the time-dependent left-moving 
fermion fields satisfy 

^-,s( w f*) = -sie~ tXl ip + - s (—v F t), (15a) 
ipl s (L + v F t) = sie- iX2 ^+.- s {L - vpi). (15b) 



Choosing t — (L + x)/v F in Eq. (p_5q), we obtain 



ip_ s (x + 2L) = sie iX2 tp + ,- s (—x). Equation (15a) gives 

■0+,s(— x) = —sie^ipl s (x) which, in combination with 

the previous equation, leads to 



^-, a (x + 2L,t) = e^_, i (s J t), 
ip +tS (x,t) = -sie lXl ipl^ s (-x,t), 



(16a) 
(16b) 



where i? = 1 + (x/ 71 ") an d X = Xi ~~ X i- Thus, we see that 



the Andreev boundary conditions (11a, lib) are equiv- 
alent to twisted periodic boundary conditions for ip- tS , 
Eq. (16a), on an interval of length twice the length of the 
original system, supplemented by the connection between 
the and "0-,s fields following from the chiral symme- 
try, Eq. ( 16b ). (Equivalently, the periodic boundary con- 
ditions can be derived for ip+,s and the chiral symmetry 
can be used to obtain ip-,s-) The problem thus becomes 
very similar to one of fermions on a ring of circumference 
2L, threaded by an effective Aharonov-Bohm flux #/2, 
the persistent current |25]-[2~7[| being the analogue of the 
Josephson current. A detailed treatment of persistent 
currents in Luttinger liquids was given in Ref. [p8| (for 
the case of spinless electrons), and we shall adopt this 
treatment in what follows. There are significant differ- 
ences between the system of electrons on a ring and an 
SLS junction, however: (i) On a ring, the number of elec- 
trons is fixed. Therefore, the system is described by the 
canonical ensemble, and is sensitive to the parity of the 
total number of fermions ]29|p8| ] . In the normal part of 
an SLS junction, the number of electrons is not fixed, so 
that the system must be described by the grand canon- 
ical ensemble, the chemical potential being maintained 
in the superconducting leads, (ii) On a ring, all four 
components of the fermion fields (^±,1,1) satisfy equiv- 
alent twisted periodic boundary conditions. In an SLS 
junction, only two of the four comp onen ts satisfy twisted 
periodic boundary conditions, Eq. (16a); once these com- 
ponents are constructed, the o ther two are found using 
the chiral symmetry, Eq. (16b). 

We now discuss the range of validity of the Andreev 
boundary conditions. The condition e <C A , which we 
needed in order to arrive at Eq. (12), means that our 



boundary conditions are capable of describing only ex- 
citations with wavelengths 1/k 3> Tiv f /Aq ~ £§■ Such 



excitations exists only in "long" junctions, i.e., L ^ £s; 
thus our treatment is valid only for this case. On the 
other hand, as follows from self-consistent calculations 
|^| , the order parameter in S gets reduced from its bulk 
value over the scale £s near the boundary, which also af- 
fects the excitations in N in the boundary region of the 
thickness of £s- Thus, the model of a step-like profile of 
A, Eq. ([!]), can adequately describe only processes tak- 
ing place in the interior of N (i.e., for x outside boundary 
layers of width £s)j where the exact shape of the pro- 
file of A in S is irrelevant. The latter condition can be 
satisfied only if L ^> £s- Therefore, the range of valid- 
ity of our boundary conditions is the same as that of 
the non-self-consistent model itself. We can also view 
Eq. ( |l2"l) as the minimal-model boundary conditions that 
describe the time-reversal process associated with An- 
dreev reflection. Thus, we now relax the assumption of 
ID superconducting leads, and regard Eqs. ( |l2| ) as the 
general boundary conditions satisfied by the low-energy 
components of the fermion fields. 



III. BOSONIZATION OF LUTTINGER LIQUID 
IN CONTACT WITH SUPERCONDUCTORS 

We now turn to the bosonization of an interacting ID 
electronic system [|30f in contact with superconductor. 
We represent the free fermion fields in the conventional 
bosonic form: 



•>P±A X ) 



1 



: exp (±iy/TT(/)± tS (x)) , 



(17) 



where a — > +0 is a convergence factor and the chiral 
bosons <j)± tS are expressed through the density (phase) 
bosons <j) s {6 s ) via 



<t>±,s( x ) = <M X ) T s (x). 



(18) 



We construct mode expansions for (f>±_ s (x) in such a way 
that the twisted boun dary conditions (16a) and the aux- 
iliary conditions (16b) are satisfied: 



3 ( x ) = ^ + ^(N s +#)£. 



<f>+A x ) 



2L 



$ a (x), (19a) 
+ 4>-„{-x). (19b) 



(The additive c-number terms have been omitted in the 
expansion for </>+, s -) Here, (p s are zero-mode operators, 
N s are operators whose eigenvalues give the winding 
numbers of the Haldane (topological) excitations 
and 4> s (x) are the nonzero- mode components of the chiral 
boson fields, which are periodic on the interval (0, 2L): 



4> s {x) 



k>0 



-ikx t 
a k,. 



Akx , 



(20) 



where k — nn/L (with n — 1,2,...), 7^. = 
exp(— an/2) / VkL, and at, s satisfy the canonical commu- 

— S ss 'Skk'- In Eqs. (19a, 191)) 



tation relations 
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the terms linear in x describe the topological excitations 
of the bosonic system, which do not conserve the total 
number of fermions. The eigenvalues of N s give the num- 
bers of fermions added to or removed from the Luttinger 
liquid. The nonzero-mode components <p s describe the 
quantum fluctuations around the topological excitations. 
These correspond to the fluctuations in the fermion den- 
sity that conserve the total number of fermions. 

We require that the chiral bosons obey the canonical 
commutation relations fl32f 

oo 

[4>±, s (x),d x '(f>±,s'(x')} = Ti$ss> ^ S ( x ~ x ' + 2nL )> 

71 — - — OO 

(21) 

where the summation over n reflects periodicity on the 
interval (0, 2L). The nonzero-mode components of, e.g., 
expansion (19a) obey 



,(a:),£W-,.s'(2')] 

OO 

= 5 ss >[ix/L + ^2 S(x- x' + 2nL)\. 



(22) 



Thus, in order for Eq. ( |2l| ) to be satisfied, the zero-mode 
operators ip s and the winding-number operators N s must 
obey 



[ip s ,N s ,]=2iS s 



(23) 



(The same result can certainly be found by consider- 
ing the commutation relations of (j)+.s-) A posteriori, 
we ca n also j ustify the choice of the coefficients 7& in 
Eqs. (19a, 19b): they were chosen in such a way that the 
commutation relations ( pl| ) are satisfied. 

Next, we introduce the charge (p) and spin (a) com- 
ponents of the boson fields: <j> p a = (0f ± </>jJ/\/2 and 
9 p / a = (6*|±0|)/\/2. The mode expansions fo r <j) u and P 
(where p = p,a) follow from the expansions 



M X ) = ■% + ^ k C0S kx { a kp + a kp 

* k>0 

<t><x{x) = \l ^M^ + i^-yk sin kx (a[ a 

* fc>0 

Op(v) = V/f ( J + % ^J J + i^2"fksinkx (al p - a 

' fe>0 

8 a (x) = "^j= + ^2 Ik cos kx (a\ a + a 
v n k>0 



(24a) 

(24b) 

(24c) 
(24d) 



where <p. p/a = ± M = (iV T - J = 

(iV| +Ni)/2, and a kp / a = (a^j ±ak,i)/v2. It is natural 
that the phase difference of the superconducting order 
parameters \i which determines the charge flow between 
the superconductors, appears only in the field associated 
with the charge current, i.e., 6 p . 



We now have to determine the topological constraints 
imposed on the winding numbers N s (and, consequently, 
on M and J). This can be don e by substituting the 
expansion, e.g., for </>_, s , Eq. (19a), into the bosonization 
formula ( p~7[ ) , and requiring that the boundary conditions 
for fermions (16a) be satisfied Pq| . [When disentangling 
the operators in the exponent of Eq. (|l7|), one must recall 
that 4> s and N s do not commute, and use Eq. (]2^).j One 
thus finds that N, satisfies 



(-) 



1, 



(25) 



i.e., that the eigenvalues of N s are odd. (Neglecting the 
operator-nature of the zero-modes and the winding num- 
bers would have led to N s being even.) Consequently, 
J+M must be odd. It is convenient to introduce an effec- 
tive win ding - number J' = J+l, so that J + d = J 1 +x/ 7r 
in Eq. (24c). Then, J' + M must be even. Comparing 
this constraint with the similar constraint on the topo- 
logical numbers in the persistent-current problem p8| , 
we see that our constraint effectively corresponds to the 
case of an odd number of fermions on the ring, in which 
case the response of the system to the twist in boundary 
conditions is diamagnetic, i.e., the free energy is minimal 
at zero twist. Tracing back through our calculations, 
we note that the diamagnetic nature of the Josephson 
current is guaranteed by the Andreev phase shift (7r/2), 
which ultimately shifts J to J + 1. The physical meaning 
of the topological constraint is quite simple: the energy 
of the LL is minimal when the left- and right-moving 
branches of the spectrum are populated symmetrically; 
changing the total number of fermions in the LL by an 
even (odd) number results in excitations with even (odd) 
total momentum quanta. 

We note that expansions similar to 



Eqs. (24a 24b , 24c , 24d ) could have been obtained by first 
deriving the boundary conditions directly for the charge 
and spin bosons from the boundary conditions for the 
fermions (in the same way that the boundary conditions 
for bosons are derived from the Dirichlet boundary con- 
ditions for fermions in Ref. |2^]), and then constructing 
expansions satisfying these boundary conditions. In this 
way, however, the zero-modes of the expansions, which 
are crucial for the topological constraints on the eigen- 
values of M and J, might have been missed. (We will 
derive and use the boundary conditions for charge/spin 
bosons in Sec. 0, when the topological structure of the 
boson fields will not be important.) 

The bosonized Hamiltonian of the LL is given by 



h \ ^ 



L/2 



dx 



K, 



■(d x <t> P ) 2 + v p K p {d x e p f 



(26) 



If the LL model originates from the Hubbard model then 
Kp/v = 1/V1 3= <?, where g = Ua/irvp, with U being 
the strength of on-site interactions and a the microscopic 
length cut-off (of order the Fermi wavelength) , and v p = 



5 



vf/K^. In addition, if the underlying SU(2) symmetry 
of the Hubbard model is intact, then K„ = 1 1331 - 



IV. JOSEPHSON CURRENT THROUGH A 
LUTTINGER LIQUID 

One of the most important consequences of induced co- 
herence in the N part of an SNS system, is the Josephson 
current through it. This current differs from its coun- 
terpart in tunnel junctions in that the critical current 
J c decays with the junction length L according to the 
power-law 1/L (for L <C Lt), rather than exponentially. 

The Josephson current in SNS junctions is affected 
strongly by the quality of the interface. The transmit- 
tance of interfaces between semiconductors and super- 
conductors varies widely, depending on the nature of the 
junction. The interface that has been studied most in- 
tensively in recent years, particularly in the context of 
mesoscopic effects, is the Nb/InAs interface. This inter- 
face is unique in the sense that a charge-accumulation 
layer is formed instead of a Schottky barrier and, as a re- 
sult, the interface transparency is quite high. (According 
to a recent measurement of the proximity effect in this 
structure, the interface transmission coefficient To f=a 0.7 
p4j.) More commonly, however, the transmittance may 
be quite low, both because of interface roughness and 
Schottky-barriers formation. Below, we calculate the 
Josephson current through the LL in t wo lim iting cases: 
perfectly transmitting interf aces ( Sec. IV A ) and poorly 
transmitting interfaces (Sec. IV B). The latter case has 
been investigated in Ref. Jl3f . 



X p/a = i^tanh 1 — 



K 2 



K 2 , 

p/o 



(29) 



We see that the phasc-diffcrcncc \ appears only in the 
topological part of 5, as it should, because the nonzero- 
mode excitations are neutral, and therefore do not con- 
tribute to the (equilibrium) charge current. We also 
note that only two of the four charge/spin bosons, viz., 
<j) a and 9 P , contribute to the topological part of £. 
(The Josephson-current problem differs in this respect 
from the persistent-current problem, in which all four 
bosons contain topological excitations.) The combina- 
tion {0 CT ,# p } commonly arises in the study of supercon- 
ductivity in LLs J33j. 

The partition function factorizes as Z = Z t (x)^n, 
where Z t / n is the contribution from the topological 
(nonzero-mode) part of £. To calculate J, we need only 
know Zt 1 which is given by 



Zt(x) = 2^ e 

J',M 



(30) 



where e p = irLrVpKp/ AvfL and e a = nv a LT / '^Ko-vfL, 
and the primed sum indicates that J' and M are con- 
nected via the constraint found in Sec. Ill (i.e., J' + M 
even). (Although the spin-part of Z t does not depend on 
X, it does not simply reduce to an overall multiplicative 
factor because of this constraint.) It is convenient to rep- 
resent the winding numbers J' and M in the following 
form: J' = 2j + kj and M = 2m + km (with kj = 0, 1 
and km = 0,1) [^8). The topological constraint is then 
satisfied for j, m = 0, ±1, . . . and kj = km- We can then 
re-write Eq. (|30|) in the unconstrained form: 



A. Perfectly transmitting interfaces 



Zt(x) = /o,p(x)/o,.(0) + /i, P (x)/i,,(0), 



(31) 



First, we consider the case of perfectly transmitting 
interfaces, in which the only scattering that takes place at 
the S /LL boundaries is Andreev reflection, single-particle 
reflection being absent. The Josephson current J is given 
by 

3 (27) 



h dx 



where 
Z is 



the 



= — fc^TlnZ is the grand potential, and 
gr and pa rtition function. Substituting 
Eqs. ( p4a| , p4bU24c1 , p4d| ) into Eq. (||) and diagonalizing 



the nonzero-mode part via a canonical transformation, 
we get for the many-body eigenenergies of the system 



£ = — 



k>0 p= 



1/2), 



(28) 



where n 



61. b,,. 



kfi = u kf ^u kfi , and the new boson operators btp, are 



connected to the old ones via a p 
in which 



bp, cosh A p 



-bl sinh \p, 



where 



/k,^(x) 



E 



1 (2n+K+ X ) 2 



(32) 



The Josephson current can now readily be calculated. 
Without writing down the exact expression (which con- 
tains, as usual for this kind of problem, Jacobi $3- 
functions |2|]), we consider only the asymptotic cases 
of low (L <g; Lt) and high (L Lt) temperatures. In 
the former case, one finds 



3 = 



eVpK p x 
L t:' 



for |x| < 7T, 



(33) 



with J{x + 27r) = J{x)- We note that the interaction- 
renormalization of the Josephson current is the same 
as that of the persistent current p8| . When the 
Luttinger-liquid Hamiltonian (|2^) is obtained as the 
long-wavelength limit of the Hubbard Hamiltonian then 
VpKp — vp. Comparing Eq. ( |33"| ) for v p K p = vp with the 
corres pon ding expressions for the non-interacting elec- 
trons gfQ, we see that the Josephson current through 



G 



the Luttinger liquid is precisely the same as through the 
non-interacting electron gas. A word of caution is neces- 
sary, however: this conclusion is only valid if backscatter- 
ing and Umklapp scattering are not taken into account. 
Even if these types of scattering are irrelevant (in the 
renormalization-group sense), they will modify the pa- 
rameters of the LL entering Eq. (Pq) , so that the equal- 
ity VpKp = v F will no longer hold |35f| . Nevertheless, the 
deviations from this equality are expected to be small. 
(For instance, in the spinless case, the maximal reduc- 
tion in the product v p K p due to Umklapp scattering is 
20%, even at half- filling, when such processes are most 
effective |^,^6|).) Also, there is a much more significant 
source of the renormalization of J which we have not 
yet taken into acc ount, viz., the non-perfectness of the 
interfaces (see Sec. IV B). 

For high temperatures (L 3> Lt) we find 



J = 8en _1 Texp (—2itaL/ Lt) 



where 



VpK p 



v F K a 



(34) 



(35) 



In the Hubbard model with the SU(2) symmetry, a — 
+ — g)/2 > 1. Thus, at high temperatures, inter- 
actions lead to the further suppression of J (in addition 
to the thermal disruption of the phase-coherence). 



B. Poorly transmitting interfaces 

Having discussed the case of perfectly transmitting in- 
terfaces, we now give a brief discussion of the case of 
poorly-transmitting interfaces. In this case, qualitative 
information can be obtained by making use of the known 
results on the interaction- induced renormalizations of the 
transmission coefficient. (For the analogous treatment of 
persistent current in imperfect LL rings, cf. Ref. J37|.) 

First, consider a non-interacting SNS "clean" system 
(i.e., the elastic mean free path being far greater than 
L), with interface transmission coefficients <C 1. 

For simplicity, we now restrict attention to the low- 
temperature case [Lt 3> L). The result for the Joseph- 
son current can be obtained from the general formula of 
Ref. p8| , Eq. (16), which expresses J through the proba- 
bility for an excitation to propagate from one interface to 
another within a certain time. Substituting the probabil- 
ity of ballistic propagation into Eq. (16) of Ref. fB8j, we 
find, after some simple algebra, that the critical current 
Jl for the structure with imperfect interfaces is 



Ja — T\TiJ c - 



(36) 



where J c is the critical current for T\ = T% = 1, which 
is given by the x~ independent factor in Eq. (|33|). In the 
interacting case, the low-transparency interfaces can be 
described within the weak-link model of Kane and Fisher 



In this model, the weak links (in our case, the in- 
terfaces) are treated as perturbations that transfer elec- 
trons between the disconnected (in zeroth order) parts 
of the LL, and the renormalization group (RG) of the 
boundary sine-Gordon model is used to find the effective 
values of the hopping (transmission) probabilities. The 
case of the double barrier has been considered in Ref. || . 
We note that: (i) because in our Hamiltonian (p6|), and 
hence in our action, the topological excitations are decou- 
pled from the nonzero-modes, the RG flow equations for 
the transmission coefficients are the same as in Ref. ||; 
(ii) the effective cut-off for the RG flow is provided in our 
case by the junction length L. Therefore, we can borrow 
the result for the renormalized product T1T2 from the 
Kane-Fisher result for the double barrier away from the 



resonance: 



T X T 2 - T&faL) 



-{l/K„+l/K a -2) 



(37) 



If the SU(2) symmetry of the underlying Hubbard model 
1, we find 



is intact, i.e., K a 



\k F Lj 



Cf a 



(38) 



which is in agreement with Ref. up to a (non- 

universal) numerical coefficient. 

Whether Eq. ([33]) or Eq. (^) is relevant to a given 
experimental situation, depends on the bare (i.e., un- 
renormalized) values of and on L. Suppose that 
Ti « Ti = Tq ps 1. Then the interface barriers can be 
treated according to the weak-barrier model [BP]. As- 
sume, for simplicity, that the potential barriers are 5- 
functions with the (bare) amplitude Vq. As Vq is small, its 
RG flow at distinct interfaces is independent, and given 
by V = Vo(L/a) il ^ K ^/ 2 . Then T is renormalized to 



T 



1 



1 



1 + (mV/hk F ) 2 



l-Tp (L\( 1 ~ K p) 



T Ka, 



For relatively short junctions, i.e., 



1-Tn 



(39) 



(40) 



the renormalization of T n due to interactions is small, and 
Eq. ( |33] ) applies. The better the interface the larger L* . 
In particular, as To L L* — > oo, in accordance with the 
previously-found virtual absence of t he re normalization 
of J for perfect interfaces (cf. Sec. IV A). For longer 
junctions, i.e., L L* , Eq. (|3q ) applies. 

Choosing g — Ua/-Kv F — 1/2 (i.e., K p « 0.8), using 
the value T = 0.7 Q and recalling that a ~ l/k F ~ 
100A in the relevant semiconductor structures, we find 
L* ~ 1 /zm. Junctions of lengths in the range 0.1 — 10 /im 
are quite common in experiments ||, so both Eqs. ( [33] ) 
and d3S| ) may be relevant in experimental situation. 

Indirectly, one also can appreciate the extent to which 
interactions renormalize the Josephson current by using 
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recent experimental results on the (dissipative) conduc- 
tance of the ultra-high mobility GaAs quantum wires 
||. As was shown in Ref. the conductance of the 
wire is reduced from the conductance quantum (i.e., e 2 /h 
per spin projection) as the temperature is lowered, the 
temperature-dependence being consistent with the the- 
ory of charge transport in dirty Luttinger liquids |}9] |IlJ . 
The absolute value of this reduction is quite small, how- 
ever: it amounts to 1 — 5% for wires of length 2 — 10 /im. 



V. 



PROXIMITY EFFECT IN LUTTINGER 
LIQUIDS 



As has been mentioned in Sec. [I A, Andreev refl.ee 



tion at the NS interface gives rise to correlations between 
electron- and hole-like excitations in N. These correla- 
tions are similar to those between the single-particle ex- 
citations in S, which can be viewed as the induction of 
superconducting off-diagonal long-range order in N due 
to the proximity of S. The presence of such order is usu- 
ally described by the (inhomogeneous) condensate wave- 
function Q F(x), defined by 



F(x) = 



(41) 



In the bulk of N, F = . The scale over which F (ex- 
ponentially) changes from its value at the NS boundary 
to zero in the bulk is given by Lt- As T — ► 0, the length 
Lt — > oo, and the exponential decay of F crosses over 
to a slower (power-law) decay. In particular, if N is a 
Fermi-liquid metal, F decays with the distance from the 
interface as 1/x (at T = 0) 0]. We now explore how 
this decay law is changed if N is in the LL state. 

Consider a semi-infinite LL occupying the half-line 
x > and connected to S at x = 0. The bosonized 
form of F(x) is given by 

F(x) = ^/e~^ e ^°)cos(VW CT (x,0))Y (42) 

where S — > +0 is a (dimensionless) cut-off parameter, 
a is the microscopic scale of the system, and 6 p (x,t) 
and (j)<j{x,T) are boson fields in the (imaginary time) 
Heisenberg representation. In Eq. (^2|), the average is 
taken with respect to Boltzmann factor e~ s l h , where 
S = S p + S a is the (Euclidean) action corresponding to 
Hamiltonian Eq. prjl), and 



S p =^ f dxdr- (d T 9 p f + v p {d x 6 p f , (43a) 



Eq. ( 11a lib ) by using the bosonized form of the fermion 
fields ([L7^ (The phase of the order parameter in S is 
now taken to be zero, as we do not consider charge-flow 
through the the interface.) Simple algebra then leads to: 



o .(O,t)=-VV4, p (O,t)=O. 



(44) 



In the (semi)-infinite geometry, the energy of topologi- 
cal excitations is infinitesimally small [jl2), and therefore 
we do not have to incorporate the winding numbers of 
such excitations in boundary conditions (p4f). As one 
might hav e anticip ated, Andreev boundary conditions for 
fermions (11a lib) impose boundary conditions only on 
those components of the boson fields that occur in the 
bosonized form of the condensate wavefunction fl42|) . 
In order to remove the divergence in Eq. (|42j ) as 

5 — > +0, we use the following trick. Consider the mod- 
ified boundary condition for the CT -field: 0<j(O,t) = 
— \/27r/4 + S. Introduce a new field cf> a = \/27r/4 — 

6 + (pa satisfying the homogeneous boundary condition 
<j) a (0,T) = 0. After this, F takes the form 



F(x) 



sin S 1 



— exp {-ir (G p (x, x, 0) + G a {x, x, 0))} , (45) 



where G p j G (x, x' ', t) is the propagator of the charge (spin) 
boson field, which satisfies 

K^(dl + v- 2 dl)G p = -S(x - x')S(r), (46) 

in which a p j a = ±1. G p obeys the following boundary 
conditions: G p (0, x',t) — 0, G p (x, x' , t)!^^^ — * and 
GJx,x',t + ft) = GJx,x',t), where (3 = 1/T. The 



Fourier transform in r of the solution of Eq. ( |46| ) is given 
by 



G p (x,x', u) 



|w| 1 sinh(|w|x<) exp(— |w|a;>), 



(47) 



where x < = min{x,x'} and x > = max{x,x'}. Inverting 
the transform, we get 



Gpix^x.O) 



1 



In (x/a) , 



(48) 



where, in order to regularize G M , we have chosen the 
same short-distance cut-off a as in Eq. (42). Substituting 
Eq. (jf§) into Eq. @, we find 



JL ( dxdT ±. {dr ^f + Va {dx< p a f . (43b ) F{x) = £ (^) 7 , with 7 = l(K a + A7 1 ), (49) 



(Note that we have deliberately expressed via those 
boson fields that enter the bosonized form of F(x).) The 
presence of S at x = imposes certain boundary condi- 
tions on these fields. We derive these boundary condi- 
tions directly from the boundary conditions for fermions 



where C is a (non-universal) numerical coefficient. In the 
absence of interactions, K p = K a = 1 and we return to 
the 1 jx scaling. In the presence of repulsive (attractive) 
interactions, 7 > 1 (7 < 1), and the condensate ampli- 
tude in the LL decays faster (slower) than in the FL. This 
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result is in accord with one's intuition: the repulsive (at- 
tractive) Coulomb interaction weakens (strengthens) the 
superconducting state induced in N by Andreev reflec- 
tion. The exponent 7 is one half of the exponent deter- 
mining the spatial decay of the (singlet) superconducting 
fluctuations in the infinite LL |33J]. 

At the first sight, the result that the profile of the 
condensate wavefunction in the LL decays faster than in 
the FL seems to contradict to the results of Sec. |v|, in 
which it was found that the junction-length dependence 
of the critical current is the same in the LL and the FL. 
Indeed, it seems natural to connect the l/x decay law 
of the condensate in the FL with the 1/L dependence 
of J c \ then, it would be reasonable to expect that the 
l/x 1 decay law of the condensate in the LL would be 
transformed into a l/L 7 -dependence of J C: even if the 
interfaces are perfect JZJ] In fact, this conclusion would 
not be valid and, as we show below, the 1/L-dependence 
of J c (at T = 0) is universal, and not connected with 
the profile of the condensate in the N region of an SNS 
junction. Consider, again, an SNS junction of length 
L » (5. Our main argument is that at T = the only 
relevant lengthscale in the problem is L; therefore, at 
distances from the interface larger than £g, the profile 
of condensate wavefunction N is described by a single 
dimensionless parameter x/L. Therefore, F(x) can be 
represented in the following form: 



F(x) = F-§-{x/L) + F+<f>+(x/L), 



(50) 



where F = F^e 1 * 1 - 2 are the values of F at x = 0(L), and 
the scaling functions <I >± (z) satisfy the following bound- 
ary conditions: $±(1) = 1(0), $±(0) = 0(1). The super- 
current flowing through the junction is given by 



J i\[F(x)^^-- C .c. 

ax 



(51) 



where A is an L-independent constant. Substituting 
Eq. @ into Eq. (§l]), we see that J can be represented 
in the form 



J = AL- l ${x/L), 



(52) 



where $(z) is a scaling function that is a combination of 
the functions $ ± (z) and their derivatives. Due to charge 
conservation, Sf{x) does not depend on x, which can be 
satisfied only if $(z) is ^-independent. Thus, we see that 
J oc 1/L regardless of the particular form of the conden- 
sate wavefunction, the latter determining only the nu- 
merical coefficient in front of the 1 /L-dependence. 



VI. SOLUTION VIA THE BOSONIZATION OF 
THE WHOLE SYSTEM 

So far, we have applied bosonization only to the 
Luttinger- liquid part of the system, i.e., to the interval 
< x < L. We can gain some further insight into our 



results by comparison with a system in which the LL oc- 
cupies the entire real line, but, by some mechanism, has 
acquired a superconducting gap when x < and x > L 
(cf. Fig. Ijc). The existence of a gap means that the usual 
Luttinger Hamiltonian is modified by the addition of the 
term 



dx\A{x)\e lx{x) + + IL 



(53) 



The corresponding (Minkowski) bosonic action is 

s = J dxdt{^ [^{d t e p f - v p (d x e p f 



+ 



2K n 



—(dt^r-v^d^y 



+\A\(:co8( X (x) + V2^(e p -<t> a )) : 

+ :cos( X (x)+V2^(e p + M) :)}, (54) 

where xi x ) is the local phase of the order parameter. In 
regions where A(x) is large, the principal effect of the 
non-linear terms is to constrain the values of p {x) and 
tp a (x) to the minima of the cosine potential, so that 



X + V2ir(6 p 
X + V2^(6 P 

Equivalently 

1 



+ 4>v) 



2mr; 
2m7r. 



/2tt 
1 

7^ 



{-X + n{n + m)); 
(+7r(n — m)). 



(55a) 
(55b) 



(56a) 
(56b) 



There are no constraints on 9 a or <fi p . The fields 9 P and <f> a 
are thus locked (modulo winding numbers) to the conden- 
sate phase in the superconducting regions. In the purely 
LL part of the system all four fields are free to fluctuate. 
The condensate therefore imposes boundary conditions 
that are essentially the same as those in Eq. (|i4|). 
The bosonized form of the number density current is 



-2v p K p 



(57) 



Substituting 8 p from Eq. (56a) into Eq. ( p7\j we find the 
T = supercurrent to be 



j(x) = 2v p K p ^-d x x- 



(58) 



We can confirm this result by considering the case of a 
Galilean-invariant system. For such a system we know 
that 



j(x) = p s v s = Ps^dxX, 



(59) 
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where p s is the density of superconducting electrons. At 
T = we will have p s = p. Comparing Eqs. ( |58| ) and 
(|59|), we see that consistency requires the equilibrium 
number density in the Galilean-invariant liquid to be 
given by 



p = 2K p v p m/Tth. 



(60) 



(The factor of 2 in this equation arises from the two spin- 
projections.) That this is correct is shown by comparing 
the commutator 



{^i,{x),—^{x')d x ^{x')] 
Ami 



-tp^(x)d x 8(x — x' 



of the charge and current in a Galilean invariant system 
with the corresponding commutator in our Luttinger sys- 
tem, viz., 



1 



[p(x),j{x')] = -2iK p v p -d x 5{x - x') 



(62) 



The Luttinger model approximates the Galilean invariant 
system by the replacement of the charge density operator 
on the right hand side of Eq. (^lj) by its expectation value. 
This confirms that E q. (6Q ) is correct. 

We now apply Eq.(|59). In the purely LL segment of 
the line (i.e., < x < L) the 8 P and cj> a fields are no 
longer constrained by the condensate. However, as we 
mentioned earlier, their values at the ends of the interval 
are fixed, just as in Eq. (Q): 



rx 2 1 
/ j(x)dx = 2v p Kp — (x2 



(63) 



This is the same result as Eq. (|35|), because the quan- 
tity found by the thermodynamic trick of differentiating 
the free energy with respect to x is the spatial average 
of the current. The advantage of Eq. (^0|) is that we 
can see that this average current is independent of the 
precise way in which the gap goes to zero as we enter 
the Luttinger link. Indeed, because the duality map be- 
tween one-dimensional charge density waves (CDW) and 
superconductors interchanges the charge- and current- 
densities, the results we have just described are just the 
dual of the well-known result in the theory of CDW sys- 
tems that the total charge induced in a region is a topo- 
logical quantity depending only on the asymptotic values 
of the CDW condensate phase fl45fl. 
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FIG. 1. (a) A Luttinger liquid (LL) conductor connect- 
ing two superconducting electrodes with phases of the or- 
der parameter xi an d X2- (J°) The model profile of the 
pair-potential used for the derivation of Andreev boundary 
conditions (Sec. |ll|). (c) Generic profile of the pair-potential 
appropriate for the bosonization of the system as a whole 
(Sec. 0). 
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